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We review the status of the calculation of the time-like splitting functions for the evo- 
lution of fragmentation functions to the next-to-next-to-leading order in perturbative 
QCD. By employing relations between space-like and time- like deep-inelastic processes, 
all quark-quark and the gluon-gluon time-like splitting functions have been obtained to 
three loops. The corresponding quantities for the quark-gluon and gluon-quark splitting 
1-^ ' at this order are presently still unknown except for their second Mellin moments. 



CLi' 1 Introduction 



The transition from coloured quarks and gluons to colourless hadrons - the so-called frag- 
mentation or hadronization process - is a Quantum Chromodynamics (QCD) phenomenon 
with many important theoretical and phcnomcnological implications for the physics at high- 
energy colliders (see, e.g., Ref. [2] for a recent summary). 

fvq I For a given fractional momentum x carried by the outgoing hadron h the fragmentation 

ff^ ■ functions (FFs) DHx,Q'^) of the final-state partons / obey evolution equations given by 
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^^ ' where the summation over j = q, q, g is understood. Q^ is a time-like hard scale such as 



the squared four- momentum of the gauge boson in the reaction e+e — >7,Z— >/i-|-X. 

The FFs are a basic ingredient at colliders, e.g., in the cross sections for particle pro- 
j^ ■ duction at high transverse momenta. However, being non-perturbative quantities, the Dl^ 

JH \ have to be determined from global analyses (fits) of various experimental hadron and photon 

production data, if possible including flavour identification. The available high-precision re- 
sults for the FFs constrain and, eventually, improve the perturbative calculation of hadron 
and photon production in a way similar to global analyses of parton distribution functions 
(PDFs) for the initial protons. With a new kinematical range accessible in proton collisions 
at LHC energies and at much increased luminosities, new opportunities for precision deter- 
minations of FFs will arise. At the same time, improved evolution equations for the FFs 
beyond the next-to- leading order (NLO) of perturbative QCD will become important. 

2 Time-like splitting functions at three loops 

The scale dependence in Eq. ^ is controlled by the time-like splitting functions P-f which 
admit an expansion in powers of the strong coupling as , 

P,n^.««(Q')) = «s/f'^(x) + alP^^^ix) + a^P^^ix) + ... , (2) 
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where we normalize the expansion parameter as Cg = as(Q^)/(47r). The three-loop contri- 
butions P-j are needed to complete predictions to next-to-next-to-leading order (NNLO) 
accuracy in perturbative QCD. 

The decomposition into flavour-singlet and non-singlet parts for the time-like evolution 
equations ^ reads 

j^^Di = Pl.^Di, (3) 



d fDs\_ ( P^. ''^l\.jDs 



th Ds EE J2{D,^+DO^ (4) 



\a/ Vqggg/ \a/ j,^]^ 

where D^^ with ^ = ±,v collectively denotes the three types of non-singlet distributions 
D^ == Dq^ ± Dq^ - {Dq^ ± Dq^) and D ^ = E"ii(^9. - Dq^). Here ® abbreviates the 
Mellin convolution written out in Eq. ([1]), while rir stands for the number of effectively 
massless quark flavors. 

Motivated by the strong similarities between the evolution equations of the time-like 
FFs in Eq. ([T]) and their space-like counter-parts for the PDFs, we have calculated in two 
publications [3,4] several of the splitting functions P^j'^ entering the perturbative expansion 
in Eq. ^ to NNLO. The leading-order (LO) terms in Eq. ^ arc identical to the space-like 
case of the initial-state PDFs, a fact often referred to as the Gribov-Lipatov relation [5] . Also 
the next-to-leading order contributions P ^ {x) have been known for long [6-8] and shown 
to be related to their space- like counter-parts by a suitable analytic continuation [6-9]. 

Following these ideas, we have derived all three-loop non-singlet functions corresponding 
to D^g and the diagonal singlet quantities Pqq and Pgg from the space-like NNLO re- 
sults computed in Refs. [10,11]. To that end, we have applied two independent methods. One 
approach relied on mass factorization, knowledge of the infrared singularities and an analytic 
continuation in the momentum fraction from x to 1/x. The other approach implemented 
the concept of universal (kinematics independent) splitting functions of Refs. [12, 13]. 

The first method started from the unrenormalized partonic structure functions in deep- 
inelastic scattering (DIS) [10,11] in dimensional rcgularization. After iteratively subtracting 
virtual contributions due to the quark or gluon form factors [14, 15] (which are well known 
also in the time-like region) we applied an analytic continuation to the x-dependent real- 
emission functions 7?.„ from x to 1/x at each order in Og similar to Refs. [6-9]. During this 
step we have also taken into account the (complex) continuation of the DIS scale q^ and the 
additional Z)-dimensional prefactor x^~^ originating from the phase space of the detected 
parton in the time-like case (see, e.g., Refs. [16,17]). The subtle point in all this is the 
treatment of logarithmic singularities for x — > 1 starting with 

ln(l — x) —>■ ln(l — .t) — ln.T + i TT . (5) 

These steps lead to the (unrenormalized) one-parton inclusive fragmentation function in 
time-like kinematics, which can be re-assembled order by order in Og, keeping the real parts 
of the continued functions TZ „ only. Finally the time-like splitting functions (and coefficient 
functions) can be extracted iteratively from the mass factorization relations. Up to two 
loops all these the steps have been checked by a direct calculation for one-parton inclusive 
electron-positron annihilation [17]. 
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At three loops we have several consistency conditions, for instance the cancellation of 
all higher poles in dimensional regularization. Most importantly, sum rules exist for the 
first moment of P-^^1 (number of fcrmions) and the second moment of -Pgg at tir = 
(momentum in pure gluodynamics). These indicate different coefRcients for the terms 

Cp Pqq (x) 71"^ In^ X of Pqq ' ' and C^ Pgg{x) TT^ In^ X of Pgg ' ' , where p^j (x) are proportional 
to the LO splitting functions. Such differences arc not unexpected in view of the imperfect 
real-virtual separation at three loops in the approach described above. 

As an independent confirmation of our analytic continuation approach including the 
sum-rule fix we thus needed a second, completely independent method. For this purpose 
we have adopted the relations between DIS and fragmentation processes of Ref. [12] which 
exploit ideas of universal, i.e., kinematics independent parton-parton splittings based on a 
physical picture of strongly ordered life-times of successive parton fluctuations. Introducing 
universal (non-singlet) splitting functions P,™'^, postulated to be identical for the time-like 
and space-like cases, Eq. ([T]) can be represented in the following manner [13] 



dlnQ2 



DL{x,Q') = / -P,rM^."^(Q')) d:J-,z'^qA , (6) 

J X ^ \ Z / 



where the notation covers both the (time- like q, a = 1) FFs and the (space-like q, Q^ = — q^, 
a = —1) PDFs. An additional shift in the argument of as in Eq. ^ also included in [13] is 
irrelevant for our discussion. 

The Os-expansion of Eq. ([6]) correctly accounts for the NLO difference S Pns (x) = 
Pns (x) — Pns (x) between the space- and time- like splitting functions [6-8]. 

At NNLO it leads to the remarkably compact expression 



SP^^^^ix) = 2 { [inx . P WJ P„(s°) + [inx . P,(°^] ® P^^^^} (7) 



where ^ = ±, v, the Mellin convolution is again denoted by ®, and 

The evaluation of Eq. ([7]) exactly coincides with the result of the above analytic contin- 
uation including the correct coefficient for the term Cp^pqq{x)n'^ In x of Pjg / , as the 
required vanishing first moment is manifest in Eq. ([7|). This result provides both the desired 
confirmation of the analytic continuation and strong evidence in support of the ansatz ^. 

The physical picture leading to Eq. ([6]) assuming strong ordering in the life-times of suc- 
cessive parton fluctuations can be applied in complete analogy to the case of gluodynamics. 
Hence a relation similar to Eq. ^ can be derived for the space- and time-like difference 
of the gluon-gluon splitting function Pgg. In particular, this reasoning again solves the 
above mentioned problem with the term (7| Pgg{x) tt^ In x of Pgg , automatically yielding 
the correct coefficient in agreement with the momentum sum rule for Us ^ 0. Moreover 
both methods, the ansatz ([6|) and the analytic continuation, agree for all remaining terms 
proportional to C|, Cjnr and CA^i (i.e., those not involving gluon emission from quarks). 

The analysis of the large- and small-x limits of the time-like splitting functions displays 
a number of interesting features. The large-x behavior of the diagonal quantities Pqq and 
Pgg is identical to that of their space-like counter-parts up to the sign of the sub-leading 
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ln(l — x) contribution as predicted in Ref. [13]. In contrast, the small-x behaviour in the 
time-hke region is markedly different as, e.g., xPJl. receives double-logarithmic contributions 
(up to a" In^" x) with large coefficients. Despite large cancellations between leading and 
subleading terms, both singlet quantities xPps and xPgg show a huge enhancement 
already at values of a; > 10"'^, sec the plots in Rcf. [4]. Although the FFs for final state 
hadrons generally do sample regions of larger x than the corresponding PDFs for initial 
protons, it will be very interesting to investigate the range of pcrturbative stability in x 
once the NNLO time-like splitting functions arc known completely. 

3 Conclusions 

We have briefly summarized the current knowledge on the time-like splitting functions in 
massless perturbativc QCD. To NLO, the complete set of functions for the evolution kernel 
in Eq. ([1]) has been known for a long time [6-8] . At NNLO accuracy, progress has been made 
by relating space-like and time-like deep-inelastic processes. Presently only the off-diagonal 
quantities Pgq and Pqg are still not completely known as a 'naive' analytic continuation 
fails, unlike at two loops, for the tt^ contributions. Access to these functions hence requires 
an extension of the methods employed so far, either by extending the concept of universal 
splitting functions [13] beyond non-singlet quantities or by refining the technique based on 
mass factorization and analytic continuation [3,4]. In particular the latter approach also 
offers other applications, e.g., to Higgs decay mediated through top-quarks as in Ref. [4] 
where highly non-trivial three-loop results of Ref. [18] were successfully confirmed in this 
manner. Thus any progress on the method of analytic continuation beyond the state of the 
art presented here may have far reaching applications. 
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